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DYNAMICAL BEHAVIOR OF A STOCHASTIC SIRS EPIDEMIC 

MODEL 

N. T. HIEU #, N. H. DU+, P. AUGER#, AND N. H. DANG# 


Abstract. In this paper we study the Kernack - MacKendrick model under tele¬ 
graph noise. The telegraph noise switches at random between two SIRS models. We 
give out conditions for the persistence of the disease and the stability of a disease 
free equilibrium. We show that the asymptotic behavior highly depends on the value 
of a threshold A which is calculated from the intensities of switching between envi¬ 
ronmental states, the total size of the population as well as the parameters of both 
SIRS systems. According to the value of A, the system can globally tend towards 
an endemic case or a disease free case. The aim of this work is also to describe 
completely the omega-limit set of all positive solutions to the model. Moreover, the 
attraction of the omega-limit set and the stationary distribution of solutions will be 
pointed out. 


1. INTRODUCTION 

The dynamics of disease spreading among a population have been investigated very 
widely in the frame of deterministic models e.g. 0,0, m, m- In such deterministic 
models, the environment is assumed to be constant. However, in most real situations, it 
is necessary to take into account random change of environmental conditions and their 
effects on the spread of the disease. For instance, the disease can be more likely to spread 
in wet (cold) condition rather than in dry (hot) condition or any other characteristics of 
the environment that may change randomly. Therefore, it is important to consider the 
disease dynamics under the impact of randomness of environmental conditions. There 
are many papers about this topic in recent years e.g. tu, HU, DU Bi- 

Weather conditions can have important effects on the triggering of epidemics. Cold 
and flu are influenced by humidity and cold temperatures. Viruses are more likely to 
survive in cold and dry conditions. Hard winds, rain, cold as well as large variations of 
temperatures are factors that weaken the immune system. Lack of sun also provokes a 
decrease of the level of D vitamin. We could also mention malaria which is influenced by 
rain and humidity level of air. Weather conditions have in general important effects on 
epidemics spread. Weather conditions change according to seasons and also to random 
variations. In the present model, we do not take into account periodic seasonal weather 
changes. Therefore, we consider an environment which is assumed to be rather constant 
on average all along the year. In this contribution, we only study the effects of random 
variations of weather conditions on the spread of epidemics. To simplify our description, 
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we also assume that only two cases can occur, bad or good weather conditions. Bad 
weather conditions corresponds to a case where the epidemics is more likely to spread 
and inversely for good conditions. Therefore, we consider that there exist two models 
associated with different parameters values corresponding respectively to bad and good 
weather conditions and that the system switches randomly from the one to the other. 
For disease models with noise, we also refer to recent contributions mm- 

The basic simplest epidemic model that we consider is the classical SIRS model intro¬ 
duced by Kernack-MacKendrick of the form (see [T5] for details) 



( 1 . 1 ) 


where the susceptible (S'), infective (/) and removed ( R ) classes are three compartments 
of the total population TV. Transitions between these compartments are denoted respec¬ 
tively by a, b , and c. They describe the course of the transmission, recovery and loss of 
immunity. 

In further studying the SIRS model, we note that the sum S +1 + R = TV and it is 
a constant of population size. So that for convenience the removed class ( R ) can always 
be eliminated. The reduction of the equation m is then 


S = -aSI + c(N — S — I) 
i = aSI — bl. 


( 1 . 2 ) 


It is easy to analyze the previous simple system El and to show that two situations 
can occur (see ng, m, E5): 

- If the basic reproduction number IZq = ^ > 1 the disease spreads among the 
population and a positive equilibrium (s*,i*) is globally asymptotically stable. It is 
therefore an endemic situation. 

- If IZo = < 1 the disease is eradicated as a disease free equilibrium (TV, 0), which 

is asymptotically stable. This situation is the eradication of the disease among the 
population. 

In this work, we shall concentrate on the switching two classical Kernack and MacK- 
endrick SIRS model, which will be chosen as the basic models for the epidemics. We 
shall assume that there are two environmental states in each of which the system evolves 
according to a deterministic differential equation and that the system switches randomly 
between these two states. Thus, we can suppose there is a telegraph noise affecting on 
the model in the form of switching between two-element set, E = {+, —}. With different 
states, the disease dynamics are different. The stochastic displacement of environmental 
conditions provokes model to change from the system in state + to the system in state 
— and vice versa. 

Several questions naturally arise. For instance, in the case where the disease spreads 
in an environmental condition, while it is vanished in the other one, what will be the 
global and asymptotic behavior of the system? Using the basic reproduction number 
1Zq of both models and the switching intensities, can we make predictions about the 
asymptotic behavior of the global system, i.e., the existence of a global endemic state or 
a disease free state? 
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The paper has 5 sections. Section 2 details the model and gives some properties of 
the boundary equations. In Section 3, dynamic behavior of the solutions is studied and 
the w-limit sets are completely described for each case. It is shown that the threshold 
A which will be given later plays an important role to determine whether the disease 
will vanish or be persistent. We also prove the existence of a stationary distribution and 
provide some of its nice properties. In Section 4, some simulation results illustrate the 
behavior of the SIRS model under telegraph noise. The conclusion presents a summary 
of the results and some perspectives of the work. The last section is the appendix where 
the proofs of some theorems are given. 


2. PRELIMINARY 


Let us consider a continuous-time Markov process , t € R+ , defined on the probabil¬ 
ity space (fl, E, P), with values in the set of two elements, say E = {+,—}■ Suppose that 

(£t) has the transition intensities + A — and —^ + with a > 0,/3 > 0. The process 
(£t) has a unique stationary distribution 


p = lim P{£ t = +} 

t->-oo 


P 

a + f3’ 


q = lim P{& 

£—>■00 


-} 


a 

a+ P 


The trajectory of (£ t ) is piecewise constant, cadlag functions. Suppose that 


0 = r 0 < Tl < r 2 < ... < t„ < ... 


are its jump times. Put 


0\ — T\ To, O'2 — Ti, ..., O n — T n T n —\... 


It is known that, if £o is given, (a n ) is a sequence of independent random variables. 
Moreover, if £o = + then a 2 n +i has the exponential density alro )00 ) exp (—at) and 
has the density /31[o,oo) exp(— Pt)- Conversely, if £o = ~ then a 2 n has the exponential 
density ad[o,oo) exp(— at) and <J 2 n+i has the density /31[o j0 o) ex P( — Pt) (see [TS] vol. 2, 
pp. 217]). Here l[o,oo) = 1 for t P 0 (= 0 for t < 0). 

In this paper, we consider the Kernack-MacKendrick model under the telegraph noise 
of the form: 

(S = -a(Z t )SI + c(Z t )(N-S-I) 

\i = a(Zt)SI-b(St)I 

where g : E = {+, —} —> R + for g = a,b, c. The noise (£t) carries out a switching between 
two deterministic systems 


and 


s = —a(+)SI + c(+)(N -S-I) 
i = a(+)SI — b(+)I, 

(2.4) 

S = —a(—)SI + c(—)(N -S-I) 
i = a(—)SI — b(—)I. 

(2.5) 


Since (£t) takes values in a two-element set E, if the solution of (12.31) satisfies equa¬ 
tion (12.41) on the interval (t„_i, T n ), then it must satisfy equation (12.51) on the interval 
(Tn) t„+i) and vice versa. Therefore, (S (r n ), I (r„)) is the switching point, that is the 
terminal point of one state and simultaneously the initial condition of the other. It is 




4 


N. T. HIEU *, N. H. DU+, P. AUGER*, AND N. H. DANG* 


known that with positive initial values, solutions to both (12.41) and (12.51) remain nonneg¬ 
ative for all f ^ 0. Thus, any solution to (22) starting in intR+ exists for all f ^ 0 and 
remain nonnegative. 


It is easily verified that the systems (12.41) and (12.51) respectively have the equilibrium 
points 


<**=>("-Ijfi) 


= n±) + c(±)' I (2 ' 6) 

and their global dynamics depend on these equilibriums. Concretely, if if > 0 then these 
positive equilibriums are asymptotically stable, i.e., when AC > lim t ^ 00 (S' ± (t),/ ± (t)) 


= (sj, i^). This is the endemic case, both susceptible and infective classes are together 
present. On the contrary, if AC ^ then lim t _ >00 (5 ,± (t),/ ± (t)) = (AC, 0) and the in¬ 
fective class will disappear. It is called the free case 


3. DYNAMICAL BEHAVIOR OF SOLUTIONS 

In this section, we introduce a threshold value A whose sign determines whether the 
system (12.31) persistent or the number of infective individuals goes to 0. Moreover, the 
asymptotic behavior of the solution is described in details. 

For any (s 0 , *o) £ intR+ with So + *o ^ A/", we denote by (5(f, so,io,u), I(t, sq, io, cu)) 
the solution of (12.31) satisfying the initial condition (5(0, so, io, w), 7(0, so,io,u)) = (so, io)- 
For the sake of simplicity, we write (5(f),/(f)) for (5(f, so, io, ui), I(t, sq, io, w)) if there 
is no confusion. A function / defined on [0, oo) is said to be ultimately bounded above 
(respectively, ultimately bounded below) by m if lim supt_>oo/M < 171 (respectively, 
liminft_ > , 00 /(f) > to). We also have the following definitions for persistence and perma¬ 
nence 

Definition 3.1. 

1) System (12.31) is said to be persistent if lim sup^^ S(t) > 0, lim sup^^ I(t) > 0 
for all solutions of (12.31) . 

2) In case there exists a positive e such that 

e ^ lim inf 5(f) ^ lim sup 5(f) ^ 1/e, 

t—yoo 

e ^ lim inf I(t) ^ lim sup I(t) ^ 1/e, 

^ t—t oo 

we call the system m is permanent. 

It is easy to see that the triangle V := {(s, i) : s ^ 0, i ^ 0; s + i ^ AC} is invariant for 
the system (12.31) . In the future, without loss of generality, suppose that ^ jjjML 

We here define the threshold value which play a key role of determining the persistence 
of the system (12.31) . 

A = p(a(+)N - &(+)) + q(a(-)N - &(-)). (3.7) 

In the first part of this section, we prove the following classification 

Theorem 3.2. 


(1) If A < 0 then linp-^oo /(f) =0 and linp-^oo 5(f) = N. 
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(2) If A > 0 the the system (12.31) is persistent. Moreover, if ^7-4 < N then the 
system is permanent. 

The proof of this theorem is divided into several propositions. Firstly, we have 

Proposition 3.3. 

a) If A > 0 then there is a 5 1 > 0 such that lim sup^^ I(t) > di. 

b) If A < 0 then lim^oo I(t) =0 and lim^oo S(t) = N. 

Proof, a) The second equation of the system (12.31) follows 

lnJ(*)-lnJ(0) = 1 _ m)d l 

Since I{t) < N, lirnsup^^ ln J(°) ^ q Therefore, 


lim sup (-j- / (a(&)N - b(£j))dt - -j- [ a(^)(N - S{t))dt) = 

t— s-oo u JO I Jo ' 

lim sup j f (a(£j)S(t) - b(^))dt ^ 0. 

t—t oo t J o 


Thus, 

lim inf -j- f (a(^ T )N - b{f£])dt < lim inf i f a{£,j)(N - S(t))dt. 
t^OO t J 0 t-HX> t J Q 

Because the process (£t) has a unique stationary distribution lim t _ >00 P{£ t = +} = p and 
lim^oo P{£ t = —} = q, then by the law of large numbers 


lim i f (a(&)N - b(^))dt = p(a(+)N - &(+)) + q(a(-)N - b(-)) = A. 

t->oo l J p 

Denote g m in = min(g(+), g{—)), g m ax = max( 5 (+),p(-)) for g = a,b,c. We have 
lim inf — [ a max (N — S(t))dt ^ lim inf — [ a(jb){N — S(t))dt 

t-s-oo t J o ' t-too t J 0 

> lim inf j f a(^){N - b(^))dt = A. 

t-s-oo t J o 


(3.8) 


On the other hand, from 

S(t) = -(a({ t )S(t)+c(t t ))I(t)+c(t t )(N-S(t)) > -(a(&)N+c(&))I(t)+c(&)(N-Sm 

it follows 


s(t) t s{0) >\f o ~(a(m+cmim+±J* c^xn- smt. 

Since lim^oo s ^~ t s ^°' ) = 0, 

limsup (j -(a(£[)N + c^ T ))I(t)dt + fj^ c (£t)( N - S{t))dtj < 0. 
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Hence, 


liminf y [ (a max TV + c max )/(t)dt ^ liminf \ [ (a(&)N + c(£,j))I(t)dt 

t-s-oo f J o t->oo f J o 

^ liminf — [ c(£A(N — S(t))dt ^ liminf i / c m i n (TV — 5(f))df. 
*->•00 t J o *->°° t J o 


(3.9) 


Combining (13.81) and (13.91) . we obtain 

1 /■* _ l 

liminf — / I(t)dt > liminf — 
t-+oo t L W 


*~>°° t J 0 ImaxN + C] 
> 


-(TV — S(t))dt 


0 u max J v ~r ’-max 
C-min 


-A > 0. 


(^max tv + 

Cmax)^n 

This inequality implies that there exists <5i > 0 such that lim sup^^ /(f) > di. 
b) From the inequality 

= a&)S(t) - b(&) < a(&)TV - &(&), 


we have 


lim sup 

i—>■ oo 


In/(f) — ln/(0) 


< lim sup i [ (a{&)S(t) - b(£j))dt < A < 0, 

t—>oo t do 


which implies that lim /(f) = 0. On the other hand, 

£—>•00 

5(f) = -a(&)5(f)/(i) + c(&)(TV - 5(f) - /(f)) > —a max NI(t) + c min (TV - 5(f) - /(f)). 
Thus, 

5(f) ^ f e “ c “ i » ( *“ t) (— 0max // + c min )/(f)df + 5(0)e _Cmint + c m i„iV f e-^^dt. 


to 


to 


We see that lim f Q e Cmin (* *)</f = . Further, by paying attention that lim /(f) = 0 

we also obtain 

r t 


£—>•00 


lim [ 

*->°° Jo 


e Cmin( * ^(-QmaxJV + C m in)/(f)df = 0. 


Hence, liminft-xx, 5(f) ^ TV. Combining 5(f) < TV for all f > 0 gets lim^oo 5(f) = TV. 
The proof is complete. □ 

By definition of A in (13.71) we have the following corollary. 

Corollary 3.4. If Ts TV then lim^oo /(f) =0 and lim^oo 5(f) = TV. 

In view of Corollary 13.41 in the following we suppose that < TV. 

Proposition 3.5. 5(f) is ultimately bounded below by 5 m i n > 0 and there is an invariant 
set for the system \2.3\) . which absorbs all positive solutions. 


Proof. Let 5 m i n be chosen such that 

- TVa(±)5 min + c(±) ( 2aU-) ~ S ' min ) > 771 > 


(3.10) 
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and let A = (S min , 0), B = (S min , N - gg, C = (gg iV - gg, D = (N, 0). In the 
interior of the triangle V we have /(f) = a(£ t )(S(t) — ^jg)/(f) < a(£ t )(5(f) — ^g)/(f) ^ 
a (&)(l^+y “ sgW) = ~ a (€t)$+jl{i) < 0 for a11 Points lying above the line BC, 
whereas S > m for all points that are below the line BC and on the left of AB by (13.101) . 
Therefore, it is easy to see that the the quadrangle ABCD is invariant under system 
(12.31) and all positive solutions ultimately go there. □ 

Corollary 3.6. If A > 0 then the system 12.31 ) is persistent. 


Proof. This result follows immediately from Propositions 13.31 and 13.51 □ 

Proposition 3.7. /(f) is ultimately bounded below by I m in > 0 < N. As a result, 

the system m is permanent. 

Proof. Since < N, we can find an 0 < £o < di such that min { — a(±)si + c(±)(AT — 

s — f) > 0 : 0 < s ^ g^y,0 < i < £o} > 0. Then, while /(f) < £o and 5(f) < gy 

T (, a(QS - b(Q)I , r , 

> —kl where k is some positive 

dl 


we have S > 0 and — = 


S —a(£ t )SI + c(£ t )(N — S — I) 


number. Denote by 7 the piece of the solution curve to the equation — = —kl starting 

dS 

at (5 m in,£o) and ending at the intersection point (giy,£i) of this solution curve with 
the line s = g^y . Let G be the subdomain of quadrangle ABCD consisting of all 
(s,i) £ ABCD lying above the curve 7 if s ^ g^y and lying above the line i = £ 1 if 


g^y < s ^ N. Obviously, G is invariant domain because 4 - > —kl, S > 0 on 7 and 

s 

i > 0 on the segment / = e 1 , < S ^ N. Since limsup/(f) > <5i > £0 and (5(f),/(f)) 

t—too 

must eventually enter the quadrangle ABCD, (5(f),/(f)) also eventually enters G which 
implies that /(f) ultimately bounded below by I m - ln = E\. □ 


Our task in the next part is to describe the w-limit sets of the system flZp . Adapted 
from the concept in [7], we define the (random) uj —limit set of the trajectory starting 
from an initial value (so,io) by 

n(s 0 ,*o,w) = nu (S(t, so, io, w), I (f, so, io, w)). 

T>0t>T 

This concept is different from the one in [5] but it is closest to that of an uj —limit set 
for a deterministic dynamical system. In the case where ff(so, io, w>) is a.s. constant, it 
is similar to the concept of weak attractor and attractor given in [2DJ HH]. Although, in 
general, the w-limit set in this sense does not have the invariant property, this concept 
is appropriate for our purpose of describing the pathwise asymptotic behavior of the 
solution with a given initial value. 

Let 7 T("(s, i) = (5 + (f, s, i),/ + (f, s, i)), (resp. irf(s,i) = (S~ (f, s, i ), I~ (f, s, i))) be the 
solution of ( 12 . dl) (resp. (12.51) 1 starting in the point (s,i) £ K+. 

From now on, let us fix an (so, * 0 ) £ R+ and suppose A > 0. This implies that at least 
one of the systems (12.41) . (12.51) has a globally asymptotically stable positive equilibrium. 
Without loss of generality, we assume the equilibrium point of the system (12.41) has this 
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property, i.e., lim^oo 7r+(s, i) = (s+,if) £ intM+ for any (s,i) £ intR+. Also, suppose 
that £o = + with probability 1. 

For e > 0 small enough, denote by Lt s (s, i ) the e-neighborhood of (s, i) and by H e C R+ 
the compact set surrounded by AB,BC , CD and the line i = e. Set 

S n {u;) = S(r n , s 0 , i 0 , w); I n (ui) = I(t„, s 0 , io, ui), Dff = a(r k ■ k < n); JFff = <j(T k -T n : k > 

We see that ( S n ,I n ) is adapted. Moreover, given £ 0j then Tq is independent of 

Tff. To depict the w-limit set, we need to obtain a key result that ( S n ,I n ) belongs 
to a suitable compact set infinitely often. In addition, we need to estimate the time 
a solution to (12.411 starting in a compact set enters a neighborhood of the equilibrium 
(s+,i+). These results are stated in the following lemmas. 

Lemma 3.8. Let J C V D {5 > 0, / > 0} be a compact set and (sf,if) G J. Then, for 
any 82 > 0, there is a T\ = > 0 such that n+(s,i) G Us 2 (sf ,if) for any t ^ T\ 

and ( s,i ) G J. 

Proof. Consider the system imp . Since (s+,*+) is asymptotically stable, we can find a 
82 = 82 ( 82 ) > 0 such that 

7r+(^ 2 (s+ *+)) cU S 2 (s+,it) Vf > 0. 

On the one hand, for (s,i) G J, Hindoo (s, i) = (s+,i+) which implies that there 
exists a T S j satisfying 

7r t + 0,i) £ Ug 2 / 2 (st,if) for all t^T si . 

By the continuous dependence of the solutions on the initial conditions, there is a neigh¬ 
borhood U S i of (s, i) such that for any (u, v) G U S i we have 

GU-g 2 (st,it). 

As a result, 

n+(u,v) G 7r+_ Ts . (Wj 2 (s+ ,i+)) C Us 2 (sf ,if) Vf ^ T si . 

Since J is compact and the family {U S i : (s,i) G J} is an open covering of J, by Heine- 
Borel lemma, there is a finite subfamily, namely { U Si u,i = 1,2, ...,n}, which covers J. 
Let T\ = maxi^j^ n {T Siii }. We see that if (s,i) G J then 7 q + (s,i) £ Us 2 (sf , *+) for any 
t^T x . □ 

Lemma 3.9. There is a compact set K G intsuch that, with probability 1, there are 
infinitely many k = k(ui) £ N satisfying (S^fc+i, hk+i) £ K. 

Proof. In the case < N, we can choose K = G, that was established in Proposition 
13.71 and see that (S^fc+i, hk+i) £ K for every k > ko. 

Suppose that Wyy < N < ^y. With (5! is shown in Proposition 13.31 by a similar 
way to the construction of the set G in the proof of Proposition ^. 71 we construct a curve 
7 + for the system (12.41) . with the initial point (5 m i n ,^i) and the end point 
then define the subdomain K of quadrangle ABCD consisting of every (s, i) G ABCD 
lying above the curve q + if s ^ ^yjy and lying above the line i = J+ in if ^yyjy ^ s ^ N. 

It is seen that K is an invariant set for the system (IO) . By Proposition 13.31 there is a 
sequence (p n ) f 00 such that 7(/z„) ^ 8 \ for all n £ N. Further, since I(t) is decreasing 
whenever S(t ) ^ ^fey, we can chose ( fj , n ) such that S(/j, n ) > ^py- If r 2k ^ /x ra < T 2 k +i, 
i.e. = +, we have (S^fc+i, l 2 fe+i) £ K because (S(p n ),T(p n )) G K and K is invariant 
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set for the system (12.41) . For T 2 k-i ^ p n < t 2 k, if SVL PT 2 k ^t<T 2 k+i ^ <^i we return to 
the above case. Otherwise, /(f) < di for any r 2 /c < f < r 2 fc +2 since /(f) is decreasing on 
[T 2 fe+i,T 2 fc+ 2 )... Continuing this process, we can either find an odd number 2m + 1 > k 
such that (5 2 m +i, / 2 m +i) £ K or see that /(f) < <$1 Vf > r 2 k- The latter contradicts to 
Proposition [373] The proof is complete. □ 


Having the above lemmas, we now in the position to describe the pathwise dynamic 
behavior of the solutions of the system (12.31) . Put 

p = j(s,i) = 71f n (n) O •• • 07T+ OTT“(s+,i+) ! 0 < fl, f 2 , • • • , f n J Tl S N| . (3.11) 

where g(k) = (-T) fe .We state the following theorem which is proved in Appendix. 

Theorem 3.10. If A > 0 then for almost all to, the closure T of T is a subset of 

n(s 0 ,io,w)- 


The following theorem provide a complete description of the w-limit set of the solution 
to (O in case A > 0. 


Theorem 3.11. Suppose A > 0, 

a) If 


i(+) b(+) c(+) 


/ \ — i/ \ (3.12) 

a(-) b(-) c(-) 

the systems dm> and m have the same equilibrium. Moreover, all positive 
solutions to the system (12.31) converge to this equilibrium with probability 1. 
b) If (13.121) is not satisfied then, with probability 1, the T = fl(so,io,cu). Moreover, 
T absorbs all positive solutions in the sense that for any initial value (so,io) £ 
intR+, the value 7 (to) = infjf > 0 : (S(t, so,*Ojw), I(t, so,io,u)) £ T Vf > f} is 
finite outside a P -null set. 


Proof, a) It is easy to see that the systems dm and (12.51) have the same equilibrium, 
(s+,i+) = ( s~,i ~) =: (s*,i*) if and only if the condition (13.121) is satisfied. Let e > 0 
be arbitrary. Since (s*,i„) is globally asymptotically stable, there is a neighborhood 
V e cW £ (s„j,), invariant under the system dm (see The Stable Manifold Theorem, [231 
pp 107]). Under the condition (13.121) . the vector fields of both systems (12.41) and (12.51) 
have the same direction at every point (s,i). As a result, V e is also invariant under the 
system (12.51) . which implies that V E is invariant under the system (12.31) . By Theorem lS.lOl 
(s*,A) £ fi(so,io,w) for almost all cj. Therefore, Ty e = inf {f > 0 : (5(f),/(f)) £ 14} < 
00 a.s. Consequently, (5(f),/(f)) £ V e Vf > Ty e . This property says that (5(f),/(f)) 
converges to (s*,?*) with probability 1 if 5(0) > 0,/(0) > 0. 

b) We will show that if there exists a to > 0 such that the point (so,*o) = 
satisfies the following condition 


det ^ + ( s ° 4 °) 4 Oo,*o)A 

\i + (s 0l io ) i~{s 0 ,i 0 )J 


± 0 , 


(3.13) 


then, with probability 1 , the closure T of T coincides H(so; w) and T absorbs all positive 

solutions. 

Indeed, let (soUo) = 7r^(s4>**~) be a point in intM^. satisfying the condition (13.131) . 
By the existence and continuous dependence on the initial values of the solutions, there 
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exist two numbers d > 0 and e > 0 such that the function p(ti,t 2 ) = 7 r t^ 7 r t 1 ( s 0 ;*o) is 
defined and continuously differentiable in (—d,d) x (—e, e). We note that 


det 

det 


dp dp 
dti ’ dt 2 J (o,o) 

-a(+)s 0 io + c(+)(N -_s 0 - *o) 
a(+)s 0 *o - 6(+)*o 


— a (—) s o*o + c(-)(A r jo - io) 
a(-)so*o - b(—)io 


^0. 


Therefore, by the Inverse Function Theorem, there exist 0 < d± < d, 0 < e\ < e such 
that p(t\,t 2 ) is a diffeomorphism between V = (0, di) x (0, ei) and U = p(V). As a 
consequence, U is an open set. Moreover, for every (s,i) £ U, there exists a (f*,^) £ 
(0,di) x (0, ei) such that (s,i) = 77 . 717 ! (so, *o) £ S. Hence, UcTc fl(so,*o,w). Thus, 
there is a stopping time 7 < 00 a.s. such that ( 5 ( 7 ), 7 ( 7 )) £ U. Since T is a forward 
invariant set and U C T, it follows that ( S(t),I(t )) £ T Vf > 7 with probability 1. The 
fact ( S(t),I(t )) £ T for all t > 7 implies that H(so,io,w) C T. By combining with 
Theorem 13.101 we get T = fl(so,*o,w) a.s. 

Finally, we need to show that when condition (13.121) does not happen, it ensures the 
existence of to satisfying (13.131) . Indeed, consider the set of all (s,i) £ intR^_ satisfying 


( S + {s,i) 
\I + {s,i) 


S 0M)\ n 


(3.14) 


or 


[— a(+)si+c(+)(N—s— «)] [a(—)s— b(— )] —[—a(—)si+c(—)(N — s—i)] [a(+)s — 6 (+)] = 0. 

(3 ' 15) 

The equation (13.151) describes a quadratic curve. However, it is easy to prove that any 
quadratic curve is not the integral curve of the system (12.51) . This means that we can 
find a to such that the point (so,*o) = 7 r^(s+,i+) satisfies the condition (13.131) . The 
proof is complete. □ 


It is well-known that the triplet (£*, S{t),I(t)) is a homogeneous Markov process with 
the state space V := E x intK^. In the rest of this section, we prove the existence of a 
stationary distribution for the process (£t, S(t) 7 1(t)). Moreover, some nice properties of 
the stationary distribution and the convergence in total variation are given. If A > 0 and 
(13.121) holds, all positive solutions converge almost surely to the equilibrium (s+,i+) = 
(s^Air)- Otherwise, we have the following theorem whose proof is given in Appendix. 

Theorem 3.12. If A > 0 and (13.121) does not hold, then (£t, S(t), I(t)) has a stationary 
distribution v*, concentrated on Ex (VDmtRt_). In addition, v* is the unique stationary 
distribution having the density f*, and for any initial distribution, the distribution of 
(£t , S(t), I(t) converges to v* in total variation. 


4. DISCUSSION 

The basic reproduction number IZo is an important concept in epidemiology. TZq is 
the threshold parameter for many epidemiological models, it informs whether the disease 
becomes extinct or whether the disease is endemic. For example, there are many recent 
papers about periodic epidemic models that concentrate on defining and computing IZo 
(see [2j , [3], [l2j, [27] [29] ). In the classic SIRS model (11.11) , 7?-o is valued by ratio - 7 ^, it 
represents the rate of increase of new infections generated by a single infectious individual 
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in a total sane population. Based on this TZq 1 we give out the key parameter A for our 
stochastic SIRS model. A reads as follows: 

A = p(a(+)N - b(+)) + q(a(—)N — £>(—)). 

This is the average of two terms associated with each system + or — weighted by the 
switching intensities. Therefore, in the stochastic model, A can be interpreted as the 
average number of infective individuals generated by a single infectious individual in a 
totally sane population for the total system with random switches. We can, therefore, 
understand that when A is positive, it signifies that asymptotically the total system will 
go towards an endemic state while the disease will vanish provided that it is negative. 
Hence, for the stochastic model, A is a very important parameter that enables us to obtain 
important informations about the asymptotic behavior of the total system. It should also 
be mentioned that in [14] , when studying the simpler SIS model with Markov switching, 
the authors also provided a similar threshold value Tq for almost sure extinction or per¬ 
sistence. In case of persistence, using the ergodicity of the Markov chain, some ultimate 
estimates for the solution were given. However, in order to study the SIRS model, we 
need to consider a system of two differential equations rather than only one differential 
equation as for the SIS model, so the method used in m may not appropriate for the 
SISR model with Markovian switching. Moreover, that method does not provide a com¬ 
plete description of the asymptotic behavior of the solution. For this reason, we treat 
our model in another way. To be more precisely, we analyzed the pathwise solution and 
obtained a better result. Moreover, some nice properties of the stationary distribution 
and the convergence of the transition probability were given. 

In this work, we have shown that the asymptotic behavior of the system depends on 
the sign of a parameter A. We have shown that when A < 0, the epidemics always vanishes 
and that when A > 0, the system is persistent leading to an endemic state. Furthermore, 
under some supplementary constraint on parameters, the system is permanent, theorem 
13.21 When A > 0, we also have shown properties of stationary distributions. We also 
illustrated our results by numerical simulations. 

We illustrate different situations in the following numerical simulations. Examples 
I and II show cases where A is positive, the first one illustrates the switching between 
two endemic systems + and —, whilst the second one depicts a system composed of an 
endemic case + and a case — for which the disease free equilibrium is stable. In both 
examples, the simulations show that asymptotically the total system persists leading to 
an endemic situation. 

The last example III considers the case of A < 0, with an endemic system + and for 
the other one — a stable disease free equilibrium. As expected, the simulation shows that 
after several switches, the disease is globally eradicated. 

Examples II and III are interesting because they illustrate a similar case, i.e. when 
systems + and — have opposite trends, system + being endemic and system — being 
disease free. In those examples, it is thus questionable to predict what will be the 
global evolution of the complete system switching at random between these two different 
situations. The answer is given by looking at the sign of parameter A which allow us to 
predict if the disease will globally invade or vanish in the long term. 

Global changes may have important consequences on the spreading of emergent dis¬ 
eases and epidemics. Therefore, it is important to provide pertinent tools that allow us 
to make suitable predictions about the possibility of emergence of a disease in a changing 
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environment undergoing climatic and environmental changes. The aim of this paper was 
to provide such efficient tools. 

As a perspective, the system would be extended to the case of a system switching 
randomly between n states, n > 2. It would also be interesting to test the model on real 
situations, like malaria, switching between wet and dry periods. Otherwise, for further 
study on the epidemic models under the effect of random environmental conditions, we 
could add some more other stochastic factors as in [5] to this SIRS model. An interesting 
perspective would also be to combine in the same model, periodic seasonal mm and 
random changes of weather conditions. 


5. Appendix 

The proof of Theorem \S.1(A With K is mentioned in Lemma ld.91 we construct a sequence 


ill = inf{2fc + 1 : {S 2 k+i, hk+i) £ K} 
r] 2 = inf{2/c + 1 > t)i : (S 2k +i, hk+i) G K} 

r)n = inf{2/c + 1 > r) n -i : (S 2 k+i, hk+i) G K} ... 


It is easy to see that {rj k =n}eJ g for any k, n. Thus the event {r)k = n} is independent 
of 5F£° if £o is given. By Lemma ITTTTT1 r] n < oo a.s. for all n. 

Let T 2 > 0,T 2 > 0. For any k G N, put Ak = {cr Vk +i < T 2 l a Vk +2 > T 2 }. We have 


P(^4fc) — P{&Vk + l < T 2, &r] k +2 > T 2 } 

OO 

= ^P{cr %+ i < T 2 ,a Vk +2 > T 2 | r]k = 2n + l}P{? 7 fc = 2n+ 1} 

n—0 
oo 

= ^ P{o"2n+2 < T 2 , cr 2n+ 3 > T 2 I T] k = 2 n + l}P{r; fc = 2n + 1} 

n—0 
oo 

= ^2 P{cr 2n+2 < T 2l a 2n+3 > T 2 }P{?^ = 2n + 1} 

n— 0 
oo 

= ^P{cr 2 < T 2 ,cr 3 > T 2 }P{^fe = 2n + 1} =P{d 2 < T 2 , (73 > T 2 } > 0. 

n— 0 
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Similarly, 


P(A fc n A k+ 1 ) — P{cr I)fc+ i < T 2 , a Vk+ 2 > T 2 , cr I)fe+1+ i < T 2 , a Vk+1+2 > T 2 } 

= P { (T %+ 1 < T 2 ,cr Vk +2 > T 2 ,a Vk+1+1 < T 2 ,cr Vk+1+2 > T 2 \ 

0^Z<n<oo 

Vk = 21 + 1, r) k +i =2 n + l}P{r; fc = 21 + 1, r; fc+ i = 2n + 1} 

= lP{ cr 2 Z +2 < T 2 ,cr 2 i + 3 > T 2 ,a 2n+2 < T 2 ,cr 2rl+ 3 > T 2 \ r) k = 21 + 1, 

0^Z<n<oo 

r/k+i = 2 n + 1 } x P{r/fc = 21 + 1 , r^+i = 2 n + 1 } 

= ^ P{d2ra+2 < T 2 ,a 2n +3 > T 2 }P{a 2 i + 2 < T 2 l a 2 i +3 > T 2 \ 

0^Z<n<oo 

r/k = 21 + 1,7y fc+ i = 2 n + 1} x P{r) k = 21 + 1 ,t^+i = 2 n + 1} 

= ^{ CT 2 < T 2 ,(?3 > T 2 }P{cr 2 i +2 < T 2 : cr 2 i + 3 > T 2 \ 

0^Z<n<oo 

r] k = 21 + 1 , r/k+i = 2 n + 1 } x P{r] = 21 + 1 , rj k +i = 2 n + 1 } 

= P{ct 2 < T 2 , tj 3 > T 2 } ^ P {°’ 2 /+2 < T 2 ,a 2 iJ r 3 > ^2 | 

0^Z<n<oo 

rj k = 21 + 1, r ] k+1 = 2n + 1} x P^*, = 2Z + 1, 7 ? fe+ i = 2n + 1 } 

OO 

= P{ct 2 < T 2 , 03 > T 2 } ^ P{ct 2;+2 < T 2 , cr 2 i + 3 > T 2 \ r] k = 21 + l}P{? 7 fc = 2Z + 1} 
i=o 

= P{cr 2 < T 2 , (J 3 > T 2 } 2 . 


Therefore, 


P(Afc U A k+ i) — 1 — (1 — P{cr 2 < T 2 ,173 > T 2 }) 2 . 


Continuing this way we obtain 


P 



1 


(1-P{<72 <T 2 ,CJ3 >T 2 }) n ~ k+1 


Hence, 

( OO OO \ 

Pi U Ai ) = : a Vn+l < T 2 , cr Vn + 2 > T 2 i.o. of n} = 1 . (5.16) 

k=± i=k ' 

Fix T 2 > 0. From j(t) = a(£ t )S(t)I(t ) - b(£ t )I(i) Js -b max I(t) and /(r,J > / min , 
it follows that /(f + r Vk ) > / m i n e _bmaxt for all Z > 0. As a result, with cr^+i < I 2 , 
T 1 1 > A *— T • P — ^ m ax^2 

Let i5 2 > 0, we choose T 2 = Ti{8 2 ) as in Lemma [3.81 for the set J = Ha- Because 
I Vk ^ <5i, it follows I Vk+1 G H a and (S Vk+2 , I Vk+ 2 ) G Wa 2 (s+,*+) provided a Vk+ 1 < 
T 2 ,a Vk+2 > T 2 . From (I5.16l) we see that {S Vk+2 ,I rik+ 2 ) G Us 2 (s+,i+) for infinitely many 
k. This means that (s+,*+) G fi(so,io 5 w) for almost all ui. 

Next, we show that ( 7 r^(s+,i+) : t ^ 0} C fZ(so,*o,w) a.s. Consider a point (s,i) = 
7 T T 3 ( s t +t)- By the continuous dependence of solutions on the initial values, for any 
J 4 > 0, there are 63^3 such that if (u,v) G Us 3 (s+ ,if) then n^(u,v) G Z^5 4 (s, i) for all 
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T 3 —T 3 <t<T 3 + T 3 . We now construct the sequence of stopping times 
Cl = inf{2fc+1 : (5 2 fe+i,/ 2 fe+i) e ^ 3 (s+,i+)}, 

C 2 = inf{ 2 fc + 1 > Ci : (S 2 k+i, hk+i) £ % 3 (s+,i+)}, 

C n = inf{ 2 fe + l > Cn-i : (S 2 k+i, hk+i) G W,s 3 (s+,*+)}... 

For (s+,i+) £ £l(so,io,u>), it follows that f n < 00 and lim ( n = 00 a.s. Since 

n—>00 

{Cfc = n} £ ?q , {Cfc} is independent of Tff. Put 

Bk = {cr Cfc+ i £ [T 3 - T 3 , T 3 + T 3 \}, k = 1,2,... 

By the same argument as above we obtain P{w : < 7 f n +i £ [ T 3 — T 3 , T 3 +T 3 ] i.o. of n} = 1. 
This implies (5<; fc +i, I^ k +i) £ Us 4 (s,i) for infinitely many times and (s,i) £ n(so,*OjW) 
a.s. Thus, { 7 r t “(s+, i+) : t ^ 0} C fi(s 0 , i 0 , w). 

Based on the continuous dependence of solutions on the initial values and using a 
similar argument, we see that { 71 ^ 0 7r t “] (s+, 7+) : t\ ^ 0,f 2 ^ 0} C f2(so,io>w). By 
induction, we conclude T C S2(so, w). Moreover, T C r2(so,*o,w) since fl(so> * 0 , w) is a 
closed set. □ 


The proof of Theorem \3.12\ We firstly point out the existence of a stationary distribution 
of the process (£ t , S(t), I(t)). From the proof of Proposition 13.31 we have 

rt 


lim inf — 

t—t OO t 


I(t)dt ^ 


(Umax A' T I'm ax )u T 


-A =: p > 0. 


Denote by 1 a the indicator function of the set A. By using the relations 

jJ 1(1)1 {/(t) < f} c ® + y ^ 


'2 ' t 

it follows, with probability 1 , that 


P N , 

H—- / 1; ru 




lim inf 

Applying Fatou lemma yields 


1 


*-«, tj 0 1 v®>v‘ a> 2N- 


Wil P { /( ^ f2AT (5 - 17) 

Consider the process (£ t , 5(f), 7(f)) on a larger state space E x (V \ {(s,i) : s = 0, 0 ^ 
i ^ iV}). it is easy to prove that (£t, 5(f), 7(f)) is a Feller process. Therefore, by using 
[22l Theorem 4.5] (or [26]) the above estimate (15.171) implies the existence of an invariant 
probability measure v for the process (£t, S (t), I (t)) on Ex (V \ {(s, i) : s = 0,0 ^ i ^ 
iV}). Since {(s,z) : * = 0, 0 < s < N} is invariant and lim t ^ oc /(t) = 0 if 5(0) = 0, 
it follows that v({(s,i) : i = 0,0 ^ s < N}) = 0. Thus, v(E x (V D intK+)) > 0. 
By virtue of the invariant property of ii'xintR.'}, the measure v* defined by v*{A) = 
u(A D E x (V D intR{_ 
v{E x (V D intM+) 

E x (V D intR{_) of the process (Ct, 5(f), 7(f)). The absolute continuity of v* is proved 


for a,nv measurable A G B(V) is a, stationarv distribution 
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as in the proof of m Proposition 3.1] while the convergence in total variation of the 
distribution of (£*, S(t), I(t)) can be referred to [10l Theorem 4.2]. □ 
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